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Abstract 

Darboux transformations in one independent variable have found numerous appli- 
cations in various field of mathematics and physics. In this paper we show that the 
extension of these transformations to two dimensions can be used to decouple systems 
of Schrodinger equations and provide explicit representation for three classes of such 
systems. We show also that there is an elegant relationship between these transforma- 
tions and analytic complex matrix functions. 
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1 Introduction 



For many years Darboux transformations in one independent variable have found numerous 
applications in various field of mathematics and physics [1-17]. (Refs. [1,2] contain an 
extensive list of references). In particular the factorization method and it generalizations 
[13-14] which have been instrumental in many physical applications (including SUSY QM 
[15] ) has its roots based on these transformations. Recently however these transformations 
were generalized and applied to systems of nonlinear equations such as the KdV hierarchy 
and soliton solutions [2,5,7]. In addition various applications of this method in geometry 
were worked out and form an important ongoing research area [2] . Extensions of the method 
to multidimensional oriented Riemann manifolds [8], time dependent potentials [9,16] and 
shape invariant potentials [17] have appeared in the literature. 

It is surprising that in spite of of this extensive research effort the theory and applications 
of these transformations in two variables and its elegant relationship to complex analytic 
function theory has not been recognized widely in the literature. [4] 

Here is a short overview of Darboux transformations for Schrodinger equation in one 
variable. (In the following we assume that all functions under consideration are smooth). 

We say that the solutions of two Schrodinger equations with different potentials u{x), v{x) 

i.e. 

0" = {u{x) + A)0 (1.1) 
iP" = {v{x) + X)ip. (1.2) 
are related by a Darboux transformation if there exist A{x), B{x) so that 



A{x)+B{x)^ 
ox 



<P{x). (1.3) 



Letting B{x) = 1 one can easily show that in order for eqs (11.11) . (11.21) to be related by the 
transformation (II. 3p A{x),u{x),v{x) must satisfy; 

A" + u + A{u-v) =0 (1.4) 

2A' + u - V = (1.5) 

Eliminating {u — v) between these equations and integration yields 

A'-A^ + u=-u (1.6) 

where u is an integration constant. Eq. (II. 6p is a Riccati equation which can be linearized 
by the transformation A = — CVC which leads to 

C" = {u{x) + v)C. (1.7) 
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Thus ( is an eigenfunction of the original eq. (11. 2p with X = u. From (11. 5p we now infer that 

v = u- 2{lnC)" (1.8) 

i.e. a Darboux transformation changes the potential function u{x) by Am = —2{ln()" where 
( is an arbitrary eigenfunction of (II. ip . 

From a more general point of view we note that Darboux transformation for equations 
in real variables form a special case of the Laplace-Darboux transformations for the complex 
variable linear equation [2] 

+ a^, + + V* = (1.9) 

(where z = x + iy). 

Our objective in this paper is to introduce a novel application of these transformations 
to decouple systems of partial differential equations in two dimensions of the form 

V^* + V(x,?/)* = (1.10) 

where = {ipi{x, y), ip2{x, y)) and V(a;, y) is a 2 x2 matrix. Such systems appear frequently 
in many applications of classical and modern physics [4,6]. 
We note that the more general form of this equation 

A{x,y)V^^ + V{x,y)^ = (1.11) 

can be transformed into (ll.lOp if A(a;, y) is (locally) invertible. 

The plan of the paper is as follows: In Sec 2 we derive the basic equations that constrain 
Darboux transformations in two dimensions. In Sees. 3, 4, 5 we solve these equations explic- 
itly for three different classes of the potential matrix V and end up in Sec. 6 with summary 
and conclusions. 



2 Darboux transformations in two Dimensions. 



We shall say that two systems of partial differential equations (PDEs) in two independent 
variables 

V2$ + U(x,?/)$ = (2.1) 

V^* + V(x,y)* = (2.2) 

are related by a Darboux transformation if there exist 2x2 nonsingular matrices with smooth 
function entries Ci(a;, y), C2{x, y), C^i^x, y) so that their solutions satisfy 



d d 
Ci{x, y) + C2(x, y)-Q^ + C3(x, y) — 



(2.3) 
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For brevity we drop in the following the dependence of the various functions on the 
independent variables. 

Using eq. (12. 3p to substitute for ^ in eq. (12. 2 p and eliminating the higher order deriva- 
tives of and using eq. (12. ip we obtain 



dC2 dCs 

dy dx 



52$ 



+ 2 



dC2 dCs 

dy dx 



dxdy 



V^Cg + 2 



dy I dy 



V^C^ + 2 ^ - C^U + VC J 



dx 



dx 



dy dx dy 



CiU + VCi ^ $ = 



(2.4) 



To satisfy this equation we treat $ and its derivatives as independent variables and let 
their coefficients be zero. This leads then to the following system of equations. 



dCo dC, 



V^Ca + 2 



V'Ca + 2 



dx dy 

dC2 ^ dCs 

dy dx 
dCi 








dx 
dCi 
dy 
dV 







C3U + VC3 = 
dU 



dC 

V^Ci - 2 — - - Cs^ - CiU + VCi = 

dy ax ay 



(2.5) 
(2.6) 
(2.7) 
(2.8) 
(2.9) 



We observe that eq. (12. 9 p can be rewritten in a symmetric form in C2, C3 in view of eq. (l2.5p . 

Eqs. (l2.5p . (l2.6p are Cauchy-Riemann equations for C2, C3. Hence the corresponding en- 
tries in this matrices must be harmonic conjugates and the entries of 

D = C2 + zC3 (2.10) 

are analytic. In view of this fact V^C2 = V^C3 = and eqs. (I2.7p . (l2.8p simplify to 

2-^-C2U + VC2 = 0, 2^-C3U + VC3 = 0. (2.11) 



dx 



dy 



From (12. lip we have 



(2.12) 
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Hence 

CsUC^^ - C3UC3 ^ (2.13) 



dx dy ^ 



Using fl2.12p to eliminate V from (12.91) leads to 

V^Ci - 2 - - Cs^ - CiU + [C3U - 2 ^]C3-^Ci = (2.14) 

ay ax ay ay 

Assuming that C2, C3 were chosen already to satisfy (12.51) - (12. 6p the system (I2.13I) -( 12TT^ 
consists of eight coupled nonlinear equations for the entries of Ci and U. However in order 
to decouple the system (12.21) by the transformation (12. 3p the resulting U must be diagonal 
or upper (lower) triangular. In the following we assume that the desired form of the matrix 
U is diagonal viz. 

jj={^^^(^^y) y (2.15) 

\ U22{x,y) J 

Under this constraint the system (I2.13p - (l2.14p is an over determined system of eight coupled 
nonlinear equations in six (unknown) functions. 

To derive solutions for these equations our strategy will be as follows: First we choose a 
proper form for the entries of the matrices C2 and C3. Then we solve (I2.13p - (l2.14p (under 
some restrictions) for Ci and U. Finally we obtain V from (I2.12p . 

In the next three sections we provide solutions and explicit examples for this decoupling 
procedure for different choices of the matrices C2 and C3 



3 C2 and C3 are Matrices with Constant Entries 

When C2 and C3 are matrices with real constants entries viz. 

C2 = ( V C3 = f ] . (3.1) 

\ d21 d22 J \ 621 622 / 

equations (12. 5p . (12.60 are satisfied by default. We also assume that Ci is of the form 

Ci= ( ^11*^^'^^ ci2{x,y) \ 2) 
\ C2i{x,y) C22{x,y) J 

In this setting eq. (12.130 leads to a coupled system of equations which can be simplified 
if we assume that 

'^21611 (^12622 „x 

621 = —; , ei2 = — (3.3) 

"11 "22 
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The resulting equations for the entries of the matrix Ci can be solved (after some algebra) 
and we find that the general solution for these entries is 



Cll = -Fii(wi), Ci2 = -^12(^2), C2l = F2i(Wi), C22 = -^22(^2), 



(3.4) 



where 

wi = diix + en?/, W2 = e22y + d22X, 

and Fij are smooth functions of the indicated variable. To proceed we assume a special form 
of the functions F^j, i,j = l,2 

h22dl2{d22X + 622?/) 



Fii = hn{dnx + en?/), Fu 
hud2iidux + eny) 



d 



22 



(3.5) 



"21 



dii 



F22 = h22{d22X + e22y)) 



where hij are constants. Using these expressions to evaluate the right hand side of 



We find: 



E 



E = V'Ci -2— ^Cg^Ci. 

dy 

V '''-'^"gr^"^^^ -2hl,{d22X + e22y)) J 



(3.6) 



(3.7) 



Substituting this result in fl2.14p and using (13. 3p we obtain only two independent equations 
for Mil and M22 



duu{x,y) duu{x,y) , ^,2 f, , ^ n 
du 7; h en h 2/1^1 (rfiix + eny) = 0, 



d' 



22" 



dx dy 
du22{x,y) , du22{x,y) 

+ 622- 



+ 2h\2{d22X + 622?/) = 0. 



dx dy 
The general solution to these equations is 

_ ^11 611(6^1 - rf^i)x^ + dl^{dl^ - el^)y'^ - 2dueu{dl^ + el^)xy 
z "11611 

^22 6^2 (ei2 - dl2)x^ + dl2{dl2 - el^)y^ - 2^22622(^^22 + e22'^)xy 



Mil 



U22 



"'22''22 



+ H2 



(3.8) 
(3.9) 



, diiy - 611X 

' dn ' 

(3.10) 
d22y - 622a;) 



d^ 



22 



(3.11) 



where Hi and H2 are arbitrary (smooth) functions of the indicated variables. 
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The general form of the matrix V (from (12.121) ) is 

/ C?11C?22(Mi1 - 2/lii) - rfi2rf2l(M22 - 2/122) C?11C?12(m22 " ^11 + 2(/lii - /l22)) \ 



V 



det{C2) 



(3.12) 



\ C?22C?2i(Mi1 - M22 + 2(/l22 - /ill)) t/2lC/l2(Mll - 2/lii) - c/iiC/22(m22 - 2/^22) / 



where det{C2) is the determinant of C2. 

Example: Let the matrices C2 and C3 be chosen as 



1 1 
1 -1 



7 '^a 



2 1 
2 -1 



(3.13) 



Furthermore let hu = /i22 = 1 and (for simplicity) Hi = H2 = 0. The potential matrix V 
for the coupled system becomes 



V 



/ -3x^ + f + 9xy + 8 -3a;2 + f + xy \ 



(3.14) 



y 3x'^ — |?/2 — a;?/ 3x2 _ |^2 _ g^^ — 8 y 

Thus the coupled system represents two coupled two dimensional oscillators. For the result- 
ing uncoupled system we have 



^ ^ \^ ~ i^^ - 5x7/ ) , M22 = -2xy. 



(3.15) 



4 The Matrix D = ^^C, where C has Constant Entries 



In this case we assume that the matrices C2 and C3 are of the form 
C, = Re(z-) I V C3 = /m(.") f 

d21 d22 J \ 621 622 



(4.1) 



where e^^dij are real and Re{z'^),Im{z^) are the real and imaginary parts of 2;", < n. 
(Similar treatment can be made if we replace by iz^). In view of f l2.5p -( l2l6l) we must have 



Since in polar coordinates Re{z^) = r" cos(n^) and Im{z"') = sm{n6) it is expedient 
to work in this coordinate system. 

It is easy to see that in this setting the right hand side of f l2.13p is zero and we obtain 
the following equation for the the elements Cij of Ci 



, ^, dcin cos((n — 1)6*) 9cj,- 
sin((n-l)^)^ ^^-^ L^^ = 0, i,j = 1,2. 



dr 



r 
7 



dd 



(4.2) 



The general solution of this equation is Cij = Fij{w) where 

cos{{n - 1)6) 



w{r, t) 



(4.3) 



and Fij are smooth functions of w. However in order to obtain a consistent and decoupled 
system of equations from fl2.14p we shall assume the following special form of Ci 



Ci 



TP ^12^22 



(4.4) 



di 



22 



With these choices for the elements of Ci the matrix E has only two independent entries 
Ell and £'22- (To derive this result we used (14. 3p ). 



— \F' + — F^l' = — 

On \- a J Hi In 



(4.5) 



where primes denote differentiation with respect to w. The matrix E takes the following 
form: 



E = — 



2n I d2ihi 



nil 

h22 



d22 



(4.6) 



di 



where ha = ha^w). Substituting these results in fl2.14p and using (14. 3 p we find that this 
system also have only two independent equations 



cos((n - 1)6)^ + r"-^ sin((n - 1)6)^ + 2nr"-i cos((n - l)e)uii = — 

or oO diiV'^"- 

If we choose Fa so that ha = then the general solution of (14. 7p is 



(4.7) 



^2n 



Hii 



sin((n - 1)9) 



(4.8) 



where Ha are smooth functions of the indicated variable. The general form of the matrix V 
in this case is 



Vii 

V12 
V21 



det{D) 
1 

det{D) 
1 

det{D) 



diid22Un - di2d2iU22 + 4r '^{d22F[-^ - 
diidi2Uii - diidi2U22 + 4r"^"((ii2F{i - 



.1 A A A 1/1 -2n/d'2ld22F[i 

«21«22Mll - a21«22M22 + 4r ( 

"11 



2n^j Tpl d2ldi2 , . 

—1 ^22J 
"22 

diidi2F22^ 
d22 

'^21-^22) 



(4.9) 

(4.10) 
(4.11) 



1 



22 



det{D) 

where det{D) = ^11^22 — (ii2'^2i- 
We present two examples. 
Example: Let n = 1 and 



d2idi2Un - dnd22U22 + 4r 



-2n 



■ d2idi2F[i 
dii 



C/11F2; 



22; 



(4.12) 



C 



1 1 
-1 1 



(4.13) 



Furthermore to satisfy fl4.5p we let Fn = F22 = -■ For these choices we find that 



A 



B 



Mil = ^, U22 = ^ 



(4.14) 



where A, B are constants. The corresponding matrix potential for the coupled system is 

A-B \ 

(4.15) 



/ A+B _ A. 

I 2r2 r6 



V 



2r2 



V 



A-B 
2r2 



A+B I ± I 
2r2 rfi / 



We observe that in this case the potentials are radial (no dependence on 6). 

Example: Let n = 2 and choose C as in the previous example. To satisfy (14. 5 p with 
ha = we choose -Fn = F22 = , q and let 



With these choices we have 



A sin 6* 5 sin 
-nil = , -"22 = • 



A sin 9 B sin 6' 

^^11 = 5 — , U22 = 



(4.16) 



(4.17) 



The corresponding matrix potential for the coupled system is 



/ (A+B) sing _ 
/ 2r' 



V 



r'^{r+cos 9)'^ 



(A-B) sin 6 
2r5 



v 



(A-B) sin < 
2r5 



(A+B) sing 
2r5 



+ 



(4.18) 



r2(r+cos6)2 



5 C2 and C3 have Orthogonal Columns 



In this case we assume that C2 and C3 have the following structure 



fii{x,y) f22{x,y) 
-fnix,y) f22ix,y) 



gii{x,y) g22{x,y) 
-9n{x,y) 922ix,y) 



(5.1) 



where fa and ga are complex conjugate. We assume also that the desired form of U is given 
by fl2.15p . Using fl3.2p to define Ci, equation f l2.13p yields after some algebra 

Mx,y)^-gn{x,y)^ = 0, A; = 1,2, (5.2) 

f22{x,y)^^ - g22{x,y)^^ = 0, k = l,2. (5.3) 

This leads us to consider the following equations 

fu{x,y)dx + gn{x,y)dy = 0, f22{x,y) dx + g22{x,y) dy = 0. (5.4) 

Although these equations are not exact an integrating factor for these equations is given by 

1 



fu{x,y) + gu{x,y) 

and 

1 



fi2{x,y) + 922{^,y) 

respectively. (This fact follows from Cauchy-Riemann equations for C2 and C3). The general 
solution of these equations can be expressed therefore by the standard formulas 

Wi(x, y) = r ^''^"^'^^ dx + r 9u{x,,y) ^^^^^^ ^^^^^ 

Jxo fiii^^y) + 9l{x,y) Jyo fii{xo,yy + 9ii{xo,yy 

It follows then that the general solution for the entries of the matrix Ci is 

cii = -Fii(wi), C12 = -^12(^2), C2i=F2i(wi), C22 = -^22(^2)- (5.6) 

However in this case we shall let -^12(^2) = F22{w2) and -^21(^1) = —Fu^wi) i.e the matrix 
Ci is of the form 

= ( ^-(^2) ) . (5.7) 

1^ -Fn(«;i) F22(u;2) ' 

The resulting matrix E is the form 

-hnix,y) h22ix,y) 



J. ^ I hii{x,y) h22{x,y) . 



where ha are defined by 

Fnjwi)" -2Fniw i)Fniw,y 
UK + 9l, 
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OfT^E) ■ 

Substituting these results in (12.141) we obtain (only) two independent equations for Ua. 

9ii^ \- lii^ \-2—-Uii = hii, 1 = 1,2. (5.11) 

ay ox oy 

Using (12.121) we find that the elements of the matrix V are 

V„^V.4(., + ..)-^^-^^. (5.12) 

V, = V., = 1k.-.„) + ^j^-^^^^^. (5.13) 

We now consider the particular case where 

^ _ . ai?e(;2") bRe{z'^) \ ^ _ ( alm{z'') blm^z"") 
^ ~ ' -aReiz"") bReiz"") J ' ^~ \ -alm{z'') blm{z'^] ' ' 

where a, b are real and n,m > 0, m ^ n. (When n = m this case reduces to the situation 
discussed in the previous section). 

As in the previous section it is expedient to work this example using polar coordinates. 
Eqs (15. 2p reduces to (14. 2 p and (15. 3p reduces to the same equation with n being replaced by 
m. The general form of the matrix Ci is given by (15. 7p where 

cos((n — 1)6*) cos((m — 1)^^) /^.^x 

^1 = ^7^31 , ^2 = . (5.15) 

Setting /ill = /i22 = we obtain from (I5.9p . (l5.10|) that 

^an(HJ±£2) 

Fii = — , or Fii = constant, (5.16) 

Cl 

tan(^^) 

F22 = — — , or F22 = constant, (5-17) 

C3 

where Cj, -i = 1 . . . 4 are constants. Eqs (15. lip for i = 1, 2 then reduce to (14.70 (with n being 
replaced by m for i = 2) and therefore (using (14. 8p ) 

«n ^ ^H. , (5.18) 
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where Hu, H22 are smooth functions of the indicated variables. The exphcit general form 
of the matrix potential V for the coupled system in this case is given by fl5.12p - fl5.13p . 

If we let n = 2, m = 1, Fa = constant, i = 1,2 and choose 



Hu 



A sin 6' 



22 



B 



(5.20) 



where A, B are constants then the matrix potential V is given by 



V 



/ A sine _|_ _B B_ _ A sine \ 



y _B _ A sing A sing _j_ ^ J 



(5.21) 



6 Summary and Conclusions. 



In this paper we showed that Darboux transformations in two dimensions can be applied to 
decouple systems of PDES. We demonstrated also the close affinity of these transformations 
to complex analytic functions. Although we were unable to solve f l2.13p -f l2.14p in general we 
were able to provide exphcit solutions for three classes of equations. Even if the decoupled 
system is not integrable analytically its numerical solution is less demanding computationally. 
Furthermore the decoupled system might provide insights about the solutions of the original 
system which are obvious directly. 
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